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Abstract. Stellar rotation produces an internal mixing of the elements due 
to shear instability and meridional circulation. This leads to observable N/C 
enhancements in massive stars above about 7-9 M© . Rotation also favours mass 
loss by stellar winds. Mass loss effects dominate for masses above 30 M©, while 
mixing dominates below that limit. The effects of mixing are also much larger at 
lower metallicity Z , because the internal O-gradients are steeper. This appears 
to be in agreement with observations in the SMC. 

At very low Z and Z = 0, mixing between the He-burning core and the 
H-burning shell leads to the production of primary N in intermediate mass stars. 
Such enrichments increase the metallicity of the rotating star, also massive Z = 
stars with moderate initial velocities currently reach break-up velocity during a 
fraction of the MS phase. Both effects favour mass loss in Z = stars, which 
have ejecta with abundance anomalies very similar to those of C-rich very metal 
poor stars. 



1. Introduction 

David Lambert has made major contributions to our knowledege of chemical 
abundances in all kinds of stars. Chemical abundances are the most decisive 
test on mixing in stellar evolution, which affects all model outputs: tracks, 
lifetimes, stellar winds, supernovae progenitors, star populations in galaxies, 
nucleosynthesis and yields, etc. Thus, the abundance determinations made by 
David Lambert are an essential piece of our astrophysical knowledge. 

Here, we examine the case of massive stars, where rotation and mass loss 
play a major role. Noticeably, several results of models at low metallicity Z also 
have implication for the abundances of very metal poor halo stars. 



2. Basic physical ingredients concerning rotation and mass loss 

We briefly summarize the phyical assumptions about the treatment of mass loss 
and rotation in stellar models. 

Structure: We cons ider the case of shellular rotation with being a function 
Sl(r) of the radius only ( Zahn 1992T ). because in a differentially rotating star the 



horizontal turbulence is strong enough to ensure a constancy of f2 in latitude. 
For differentially rotating star s, the structure equations need to be properly writ- 
ten (|Mevnet fc Maeder 1 997 h which is often not the case in current literature. 



Structural effects due to the centrifugal force are in general small in the interior, 
while the distorsion of stel lar surface may be large enoug h to produce significant 



shifts in the HR diagram (|Maeder Pevtremann 197 
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Internal transports of chemical elements and angular momentum: For these 
two transports, we consider the effects of shear mixing, of meridional circulation 
and their interactions with horizontal turbulence. The equation of evolution of 
chemical abundances due to the trans port and nuclear processes is at each level 
in lagrangian coordinates ( Zahn 1992? ) . 
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The equations of transport of angular momentum in the vertical direction is 

4 (*«)«, - hi + hi ^ ) ■ < 2 > 

f2(r) is the mean angular velocity at level r. U(r) is the radial term of the 
vertical component of the velocity of the meridional circulation. -D s hear is the 
coefficient of shear diffusion. From the two terms in the second member in this 
last equation, we see that advection and diffusion are not the same. Contrarily 
to what is often made in recent literature, the transport of angular momentum 
by circulation cannot be treated as a diffusion, (by doing so one may even have 
the wrong sign for the transport of angular momentum !). For the changes of the 
chemical elements due to transport, we see that the diffusion coefficient is the 

sum of -Dshear and Deg = ^g^p • D e g expresses the resulting effect of meridional 
circulation and of a large horizontal turbulence ( Chabover Sz Zahn 19921 1. This 



expression of D e g tells us that the vertical advection of chemical elements is 
inhibited by the strong horizontal turbulence characterized by D\. The usual 
estimate of the coefficient of horizontal turbulence was = —r \2V(r) — 
aU(r)\ (IZahn 1992ft . More recent estimates gi ve higher values for this coefficient 
( Maeder 2 003). confirmed by further studies (|Mathis et al. 20041 ) . The diffusion 



coefficient -D s hear by shear mixing essentially behaves like 
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with p=0.884. Dshear ls also modified by the horizontal turbulence D^, in the 
sense that a larger leads to a decr ease of the mixing. Th e radial component 
U(r) of the meridional circulation is ( Maeder Sz Zahn 1998ft . 
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is the reduced mass, with the notations given in the 



quoted paper. The driving term in the square brackets in the second member 



is Eq. It behaves mainly like Eq 
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average on the considered equipotential. The term with the minus sign in the 
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Figure 1. HR diagram of massive stars with Z=0.02, for rotating stars with 
^ini = 300 km/s (continuous lines) and for non-rotating stars (dotted lines). 



square bracket is the Gratton-Opik term, which becomes important in the outer 
layers when the local density is small. This term produces negative values of 
U(r), meaning that the circulation is going down along the polar axis and up in 
the equatorial plane. This makes an outward transport of angular momentum, 
while a positive U(r) gives an inward transport. At lower Z, the Gratton-Opik 
term is negligible, which contributes to make larger O-gradients in lower Z stars. 

Mass loss in rotating stars: Rotatio n certainly influences the mass loss rates 
by stellar winds (jHeger fc Langer 19981 '! . For a rotating star, one must consider 
the flux F(J)) at a given colatitude as given by von Zeipel's theorem. Thus, 
in a rotating star, the Eddington factor becomes a local quantity T^-i?). We 
define it as the ratio of the local flux F($) given by the von Zeipel theorem to 
the maximum possible local flux, which is -Fii m ($) = — ^py5cff ($)• Thus, one 
has (jMaeder fc Mevnet 2000h 
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where the opacity k("&) depends on the colatitude since T e g- also depends on 
'&. This shows that the maximum luminosity of a rotating star is reduced by 
rotation. It is to be stressed that if the limit ro(i?) = 1 happens to be met in 
general at the equator, it is not because <? e g is the lowest there, but because the 
opacity is the highest! 

Often, the critical velocity in a rotating star is written as v^ rit = ^^(1 — T). 
This expression is incorrect, since it would apply only to uniformly bright stars. 
The critical velocity of a rotating star is given by the zero of the equation 
expressing the total gravity qw. = <feg + q r lr\ = 5grav + 9mt + 5md- This equation 
has two roots (jMaeder Sz Mevnet 20001 ). The first that is met determines the 
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Figure 2. Left: Evolution during the MS phase of the angular ve locity 
in the interior of a 20 M star at Z=0.02 ijMevnet fc Maeder 200(3) . Right: 
the same for a 20 Mp, star at Z = 10~ 5 . T he values of the central H-content 
X c are indicated ijMevnet fc Maeder 200l . 



i 



critical velocity. The first root is as usual f cr it,i = (§7f~) 2 > where R p b is the 
polar radius at break-up. The second root t> cr it,2 applies to Eddington factors 
bigger than 0.639. It is equal to 0.85, 0.69, 0.48, 0.35, 0.22, times u critl i for 
r = 0.70, 0.80, 0.90, 0.95, 0.98 and 1.00 respectively. 

The theory of radiative winds applied to a rotating star leads to an expres- 
sion of the mass flux as a function of colatitude. We may distinguish two effects, 
a) The "g e fr~effect" : the higher effective gravity makes a higher T e g- at the pole 
and favours polar ejection. For a star hot enough to have electron scattering 
opacity as the dominant opacity source from pole to equator, the iso-mass loss 
curve has a peanut-like shape, b) The "opacity-effect" : if the T e g of the star is 
lower than about 25 000 K, a bistability limit i.e. a steep increase of the opacity 



( Lamers et al. 19951 ) may occur somewhere between the pole and the equator, 



due to the decrase of T e fj from pole to equator. This "opacity-effect" produces 
an equatorial enhancement of the mass loss. The anisotropies of mass loss influ- 
ence the loss of angular momentum, in particular polar mass loss removes mass 
but relatively little angular momentum. This may strongly influence the further 
evolution. We may estimate the mass loss rates of a rotating star compared to 
th at of a non-rotating star at the same location in the HR diagram. The result 
is (jMaeder fc Mevnet 20o6T l 
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where T is the electron scattering opacity for a non-rotating star with the same 
mass and luminosity, a is a force multiplier ( Lamers et al. 19 95). For B-stars 



far from the Eddington Limit, —1.5 For stars close to T = 1 the increase 

of the mass loss rates may reach orders of magnitude. 
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Figure 3. 3-D diagram showing the HR diagram horizontally and N/H on 
the vertical axis. The values of N/H in number are to be multiplied by 10, 
they are normalised to the value on the ZAMS. Case of a 15 M model with 
t>ini = 300 km/s. 



On the whole, there are 3 cases of stellar break-up: 1.- The r-Limit, when 
radiation effects largely dominate; 2.- The f2-Limit, when rotation effects are 
essentially determining break-up; 3.- The fiT-Limit, when both rotation and 
radiation are important for the critical velocity. 



General results 



Grids of models have been made (|Mevnet fc Maeder 200lh at Z = 0.020, 0.004 



and 10 . Fig. 1 illustrates the HR diagram for non-rotating and rotating stars 
at solar metallicity. The rotating models have an initial velocity Vi Q i of 300 km/s, 
which gives an average velocity of 220 km/s during the MS phase, corresponding 
to observations. We notice the following points: 

• Rotation increases the MS lifetime with respect to non-rotating models 
(up to about 40 %). 

• The values assigned from isochrones with an average rotation velocity typ- 
ically lead to ages 25% larger than without rotation. 

• Rotation strongly affects the lifetimes as blue and red supergiants (RSG). 
In particular in the SMC, the hig h observed number of RS G can only be 



accounted for by rotating models (M evnet fc Maeder 2 001). 
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Figure 4. Left: Evolution of the ratio N/C (in number) of nitrogen to 
carbon for 9 M models of different Z and velocities. The ratio are normalized 
to their initial values. The long-dashed line at the bottom corresponds to a 
non-rotating model of 9 M Q at Z = 10~ 5 . Right: the same at the surface of 
9 and 40 M Q models of different Z. 



• Steeper gradients of internal rotation Q are built at lower Z. Fig. 2 shows 
the evolution during the MS phase of Q in models at Z = 0.02 (left) 
and at Z = 10~ 5 (right). The steeper ^-gradient at lower Z favours 
mixing. There are 2 reasons for the steeper fi-gradients. One is the higher 
compactness of the star at lower Z. The second one is more subtle. At 
lower Z, the density of the outer layers is higher, thus the Gratton-Opik 
term is less important. This produces less outward transport of angular 
momentum and contributes to steepen the f2-gradient. 

• At lower Z, rotating stars more easily reach break-up velocities and may 
stay at break-up for a substantial fraction of the MS phase. 

• There are various filiation sequences for massive stars, as shown below. 

M > 9OM : O - Of - WNL - (WNE) - WCL - WCE - SN (Hypernova low Z ?) 
60 - 90 M : O - Of/WNL^LBV - WNL(H poor)- WCL-E - SN(SNIIn?) 
40 - 60 M : O - BSG - LBV & WNL -(WNE) - WCL-E - SN(SNIb) 

- WCL-E - WO SN (SNIc) 
3O-4OM : O - BSG - RSG - WNE - WCE - SN(SNIb) 

OH/IR <=> LBV ? 

25 - 30 M : O -(BSG)- RSG - BSG (blue loop) - RSG - SN(SNIIb, SNIIL) 
1O-25M : O - RSG - (Cepheid loop, M < 15 M ) RSG - SN (SNIIL, SNIIp) 

The sign <5 means back and forth motions between the two stages. The limits 
between the various scenarios depend on metallicity Z and rotation. The various 
types of supernovae are tentatively indicated. 
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Figure 5. Left: Evolution as a function of the actual mass of the abun- 
dances (in mass fraction) at the surface of a non-rotating (upper panel) and 
a rotating (lower panel) 60 M Q stellar model. Right: Evolutionary tracks 
in the X s versus log L/ Lq plane, where X s is the hydrogen mass fraction 
at the surface. The initial masses are indicated. Long-short dashed curves 
show the evolution of Z = 0.004 models, dashed-dotted curves, short-dashed 
curves and continuous lines show the evolutions for Z = 0.008, 0.020 and 
0.040 respectively. 



4. Evolution of surface abundances in massive stars 

4.1. N/C enrichments during the MS phase at various metallicities 

As a result of the rotational mixing, products of the CNO processing are reach- 
ing the stell ar surface during MS evolution and produce N/C enhancem ents, 
as observed ( Gies fe Lambert 19921 : iHerrero et al. 19921 : iLvubimkov 1996( ) since 



long in OB-stars. Fig. 3 illustrates in a 3-D plot the growth of N/H on the 
track of a 15 Mq rotating with an initial rotation Wi n ; = 300 km/s. This velocity 
corresponds to the average observed velocity during the MS phase of OB stars 
of about 220 km/s. We see that the N-enrichment progressively occurs during 
the MS phase, then it keeps about constant during the crossing of the HR and 
rises up again due to the convective dredge-up in red supergiants. 

Fig. 4 left shows the evolution of the NJ C ratios in models of rotating stars 
with 9 Mq for initial Z = 0.02, 0.004 and 10~ 5 . At zero rotation, there is no 
enrichment during the MS phase (except at .2=0.02 for M > 60 Mq due to very 
high mass loss). At 9 M for solar Z and an initial rotation of 300 km s _1 , 
the N/C ratio increases by about 0.4 dex during the MS phase. The relative 
values of the N/C ratios increase with decreasing Z, in particular N/C increase 
by two orders of a magnitude for Z = 10~ 5 . This results from the steeper 
fi-gradients and greater compactness. Of course, large N/C enhancements are 
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Figure 6. Variations of the number ratios of SN types lb and Ic with re- 
spect to SNII. The crosses with error bars correspond to observational data 
l|Prantzos fc Boissier 2003ft . The dotted line is a fit proposed by these au- 
thors. The continuous and dashed lines show the predictions of rotating and 
non-rotating models. 



accompanied by small enrichments in helium, typically of a few hundredths. The 
larger N/C enrichments at lower Z have been nicel y confirmed by abundances 
determinations in the SMC ( Venn &; Przvbilla 2004T ) . They found relative N/ H 



enrichments for A-type supergiants in the SMC reaching an order of magnitude 
or more, while in the Milky Way the enrichments are only by a factor 2 to 3. 

During the He-burning phase, the rotation velocities become all the same 
whatever the initial rotational velocities. Thus, in the He-burning phase we may 
have very different surface chemical compositions for actually similar rotation 
velocities. 

Fig. 4 right compares the relative N/C enrichments for initial 40 and 9 Mq 
stars. For high mass stars (> 30 Mq), the higher Z models have the larger 
enrichments, i.e. a situation opposite to that described above; the reason is that 
mass loss effects dominate over mixing. Higher Z stars have a higher mass loss, 
which peels off the stars and make the products of CNO process visible at the 
stellar surface. 

4.2. The abundances in Wolf Rayet stars 

For WR stars, the first major constraint is to explain their relative number ratios 
to O-type stars, which strongly decrease for lower Z. Rotation makes massive 
stars to enter earlier in the WR phase and this also increases the lifetimes as 
WR stars. When account is given to both rotation and to the dependence of the 
mass loss rates to metallicity Z, th e comparisons with observations are greatly 
improved ( Mevnet Maeder 2004h . 



Fig. 5 left shows the evolution of surface abundances for a 60 Mq star with 
and without rotation. Rotation makes smoother changes of abundances, due 
to internal mixing. For WN stars, the transition phase, when they still have H 
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Figure 7. Evolutionary tracks for models at solar and Z = l|Fei)oo 1999th 



present, becomes longer due to rotation and this increases the late WN phase 
(WNL), where H is usually present. The CNO abundances at the end of the 
WNL phase are the same for rotating and non-rotating models, because they are 
model independent and determined just by CNO nuclear equilibrium. The same 
is true for the early WN (WNE), which in general have no or little H present. 
Indeed, CNO abundances in WN stars provide a unique test of the physics of 
the CNO cycle. 

For the WC phase, the milder composition gradients, when revealed at the 
surface, make smo other transitions with lower C/He ratios, in good agreement 
with observations ( Crowther et al. 19951 ). The abundances in WC stars are not 
equilibrium values, but are products of the partial He-burning, thus they are 
model dependent and offer a most interesting test on the stellar models. Thus, 
the abundances in WN and WC stars tell us different stories. 

Fig. 5 right shows the evolution of the H-surface content X s vs. luminosity. 
This is a very constraining diagram especially for the transition stages from 
O, Of, LBV to WN stages. It is useful for establishing the proper filiations 
between such stars. It clearly supports the view that there may be back and 
forth evolution between LBV and WNL stars, and that in general WN stars 
succeeds the LBV stage. We also see that descendants from high masses at 
higher Z significantly decrease in luminosity during their evolution. 

The rare WO stars, characterized by a high O/C ratio represent a more 
advanced stage of nuclear processing. Curiously enough, such stars which may 
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Figure 8. Variation as a function of the initial mass of the stellar yields in 
14 N for different metallicities and rotational velocities. The continuous lines 
refer to the mo dels at Z = 10~ 5 , the dot ted lines show the yields from models 
at Z = 0.004 ijMevnet fc Maeder 2001|) . the dashed lines present the yields 
for solar Z models. The filled squares and circles indicate the cases without 
and with rotatio n respectively (v\r,\ = 300 km s^ 1 ). The crosses are models 
at Z = 0.1 Zr* llWooslev fc Weaver 1995^ . the empty squares at Z = 0.004 
llvan den Hoek fc Groeneweeen 1997f) and th e empty triangles are for solar Z 
models l| van den Hoek fc Groeneweeen 1 997) up to 8 Mq with complements 
<|Wooslev fc Weaver 1995^ 



be the progenitors of supernovae SNIb/c are found only at lowe r Z. The physi- 
cal reasons of that have been explained ((Smith &i Maeder 199 lh : lower Z stars 
become WC stars (if they do it) only very late in evolution, i.e. with a high 
O/C ratio. At the opposite, at higher Z the WC stars may occur at an early 
stage of He-processing, i.e. with a low O/C ratio. Fig. 6 comapres the observed 
and predicted variations of the number ratio of SNIb or c and SNII and shows 
an interesting agreement. This is noticeable because of the possible connection 
WO stars - SNIb/c - GRBs (7 Ray Bursts). In this context, we consider WO 
stars as good candidate for GRB progenitors. 



5. Abundances and yields in very low Z and Z = stars 

Evolutionary models at Z = behave very differently from models at solar 
composition, as shown by Fig. 7, in particular th e evolutionary scenarios and 
final stages are different ()Heger &: Wooslev 2 002). Let us concentrate on the 
chemical abundances. For metallicities lower than that of the SMC, the internal 
mixing in intermediate and massive stars is sufficient to bring new C from the 
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Figure 9. Chemical composition of the matter ejected through stellar winds 
by a 60 M© at Z = 0.00001 with an initial velocity equal to 2/3 of the break- 
up velocity. Note that log (Zq/Z)=3.3 in this diagram. The horizontal dotted 
lines show the solar ratios (middle) and a factor two enhancement/depletion. 



He-burning core to the H-shell where it is then turned to N. This nitrogen 
is called "primary" as it does not result from CNO elements initially present. 
Fig. 8 shows the yields in N as a function of the initial mass for Z = 0.02, 0.004 
and 10" 5 . At Z = 0.02 and Z = 0.004 the level of N production is just that 
resulting from the initial CNO elements, even if at Z = 0.004 the N-enrichment 
is somehow larger than at Z = 0.02 (cf. Sect. 4.1). However, at Z = 10 -5 
there is an overproduction of N by 1 or 2 orders of magnitude. This primary N 
is mainly produced in intermediate mass stars, but the contribution of massive 
stars is also significant. 

Remarkably, at very low Z the fast rotating stars of intermediate masses 
which reach the TP-AGB phase (this occurs for M < 7 M©) get a higher Z 
during this phase due to their enrichment in CNO elements. As an example, a 7 
M has X(CNO) = 3.1 TO -3 which is 430 times the initial CNO content. Thus 
very low Z stars may get a higher Z value in post MS phase, which may drive 
large mass loss rates. 

Intermediate and massive stars of low Z easily reach break-up velocities 
during a significant fraction of their MS phase. This results from the growth 
of the rotation velocity at their surface, due to the very weak mass loss and to 
the mild internal coupling (due to meridional circulation) with the contracting 
core, which spins faster and faster during evolution. As a result of these 2 effects 
(mass losses at break-up and from higher surface Z), massive rotating stars of 
very low Z also lose a large fraction of their masses. As an example, a star 
with an initial mass of 60 Mq at Z = 10~ 5 finishes its life as a 30 Mq core 
and may possibly lead to an hypernova forming a black hole. The composition 
of the wind ejecta of this 60 Mq model is illustrated in Fig. 9. It shows large 
overabundances of C, N and O, as well as lower enrichments in Ne, Na and 
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Figure 10. Relation between the final and the initial mass for rotating mod- 
els at various Z. The final masses of models at very low and zero Z are still 
uncertain, but they probably significantly differ from the line of those with 
zero mass loss rate M = 0. 



Al from secondary cycles associated to the CNO cycles. Remarkably, these 
abundance anomalies well correspond to those observed at the surface of the 
very metal poor halo stars showing C excess. If there is some other ejecta at the 
time of the supernova explosion, this material with anomalous abundances will 
be diluted in the current ejecta with a-rich nuclei. 

Thus, we see that if very low Z stars rotate at a significant rate, they may 
interestingly contribute to account for some abundance anomalies observed in 
very metal poor halo stars. In this respect, it is interesting to mention that the 
the ratio of Be-stars to a ll B stars is 3 to 4 times higher in the SMC than in 
the solar neighbourhood (M aeder et al. 1 999). This suggests that the fraction 
of fast rotating stars is higher at lower Z. 



6. Chemical yields from rotating stars 

The relations for various metallicities Z between the final masses at the time 
of supernovae explosion and the initial stellar masses are still uncertain due to 
the uncertainties in mass loss rates. The present mass loss estimates lead to 
the relations illustrated in Fig. 10 for Z > 0.004. The models are based on the 
expression of the mass loss rates as a function of metallicity given by M{Z) = 
M{Z Q ) {Z/Z Q ) a (|Kudritzki 2 002b a relation which is certainly different in the 
various evolutionary phases, since the mechanisms of mass loss are not the same. 
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Figure 11. Left: The chemical yields for models w ithout rotation. R ight: 
The chemical yields for models with Vi ni — 300 km/s (Hi rschi et al. 2004|) . 







OJ 




a e 




"3 
m 




"ST 




s 






Total\, 


X 4 


(SN + Wind) 


VI 




-a 




yie] 


4 He 


lar 




Lcl 




m 











) 10 


Fi 


gure 12 



km s~ 



16 







40 &0 BO 



Initial Mass [M ] 



ir Total 




50 60 70 80 90100 



Initial Mass [M ] 



For the very low and zero Z, the uncertainties are still too large, so that we do 
not know the final masses. However, we emphasized above that, in very low 
Z stars, surface enrichments due to internal mixing may enhance the mass loss 
by stellar winds, also the stars may reach break-up and experience significant 
mass loss. The values of the final masses are critical for determining the types of 
supernovae and also for nucleosynthesis, because as is evident what is escaping 
in the winds is not further nuclearly processed. 

The mode l evolution of rotat ing stars has been pursued up to the presu- 
pernova stage dHirschi et al. 2004T). since we know that nucleosynthesis is also 
influenced by rotation ( Heger et al. 200(1 ) . Fig. 11 show the chemical yields 
from models without and with rotation. Fig. 12 shows these yields multiplied 
by the initial mass function (IMF). The main conclusion is that below an initial 
mass of 30 M@, the cores are larger and thus the production of a-elements is 



14 



Maeder et al. 



enhanced. Above 30 M , mass loss is the dominant effect and more He is ejected 
before being further processed, while the size of the core is only slightly reduced. 
When we account for the weighting by the IMF, the production of oxygen and 
of a-elements is globally enhanced as illustrated by Fig. 12, while the effect on 
the He-production in massive stars remains limited. It will be interesting to 
explore the consequences of these new yields for the chemical evolution of the 
Galaxy. 

As a general conclusion, we see that the chemical abundances at all stages 
are a most constraining test of stellar evolution and in particular of the internal 
mixing process and of mass loss, which are effects strongly influencing all the 
outputs of evolution and nucleosynthesis. 
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